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Abstract 

We obtain sharp U' and Holder estimates for the Neumann function of the operator 
V • 7V — ifc on a bounded domain. We also obtain quantitative description of its 
singularity. We then apply these estimates to quantitative photo-acoustic imaging in 
inhomogeneous media. The problem is to reconstruct the optical absorption coefficient 
of a diametrically small anomaly from the absorbed energy density. 
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1 Introduction and motivation 

The purpose of this paper is to derive sharp estimates of the Neumann function of the 
operator V • 7V — ik and its derivatives, where 7 is an (scalar) elliptic coefficient defined on 
a bounded domain 51 C R'' > 3) and A: is a positive constant. The Neumann function of 
V • 7 V — ik inQ, is the function iV:51xri^CU{oo} satisfying 

{-{^■^^-ik)N{;y)=5y in 

y -/VN{-,y) -n^O on dn, ^ ' ' 

for all 2/ S ri, where Sy is the Dirac mass at y and n is the outward unit normal vector 
field on 9r2 (see subsection 12.21 for a precise definition of the Neumann function). The 
function N{x,y) has a singularity at x = y. We are particularly interested in describing in 
a quantitative manner the singularity of N{x,y) and its dependence on the parameter k. 

The investigation of this paper is motivated by quantitative photo-acoustic imaging, 
particularly by the recent work [5]. 
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The purpose of quantitative photo-acoustic imaging is to image the optical absorption 
coefBcient from the absorbed energy. The absorbed energy is obtained from boundary 
measurements of the pressure wave induced by the photoacoustic effect. We refer to [Tj and 
references therein for recent development on this inverse problem. Reconstruction of the 
optical absorption coefficient, /Iq, from the absorbed energy, A, is more delicate than the 
reconstruction of the absorbed energy from the pressure wave since fia is related to A in 
an implicit and no n- linear way (see Section |3]). One direction of research in quantitative 
photo-acoustic imaging is to reconstruct the absorption coefficient of diametrically small 
unknown anomalies. In [51 [3], efficient methods to reconstruct fia from A are proposed 
and implemented numerically when there is a small absorbing anomaly in the background 
medium. The methods use in an essential way an asymptotic expansion of A in terms of 
when the diameter of the anomaly tends to 0. The asymptotic expansion is derived using 
estimates of the Neumann function under the assumption that the scattering coefficient of 
the medium is constant. In order to extend the results of [21 [3] to inhomogeneous media, 
we shall derive sharp estimates of the Neumann function of problem (jl.ip . which is exactly 
what this paper aims at. 

To describe the kinds of results obtained in this paper, let us fix a point z € ft {z indicates 
the location of the anomaly), and let 7* := j{z). Let r(a;) :— — l/(47r|x|) be a fundamental 
solution of the Laplacian in three dimensions. Then, we will show by precise estimates 
depending on k that the singularity of N{x, z) for x near z is of the form -^T{x — z). We 
also show that the singularity of the derivatives of N(x, z) is given by the derivatives of 
-^r{x — z). We also derive L^, pointwise, and Holder estimates of the Neumann function 
iV. We then use these estimates to derive an asymptotic expansion in inhomogeneous media 
where the scattering coefficient /is is not constant. 

This paper is organized as follows. In Section 2, we derive and pointwise estimates 
of the Neumann function N. In Section 3, we show how these estimates can be used for 
reconstructing the absorption coefficient of a small absorbing anomaly. 



2 Estimates for Neumann functions 

This section is devoted to the study of the Neumann function for the operator L given by 

Lu = V • (7VU) — iku 

in a bounded domain il C M** with > 3. Here, we assume that fc is a positive constant 
satisfying k > ko for some ko > and 7 : il — > M satisfies the uniform ellipticity condition 

1^ < -/{x) < yxen, (2.1) 

for some constant f G (0, 1]. 

We first introduce some (standard) notation and definitions that will be used throughout 
the paper. Let G M'* (d > 3) be a bounded Lipschitz domain. We call diam(n) the least 
upper bound of the distances between pairs of points in ft. We say that a function / on Q, 
admits a modulus of continuity 9 if 6 : — >■ R-)_ is a nondecreasing function such that 

\.fix)-f{y)\<ei\x-y\), Vx,yen. 

For < A < 1 and / G C'^^^(il), we let [/]o.A:f2 denote the A-H61der seminorm of / in il; i.e., 

\f{x)-f{y)\ 
[Jlo.x-.n = sup — j — . 

x,yen:x^y F — y\ 
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For p > 1 and m a non- negative integer, we define the space W^"''P(ri) as the family of 
all m times weakly differentiable functions in LP(r2), whose weak derivatives of orders up to 
m are functions in Lv{n). We let W^'^^iVl) to be the closure oiC^{n) in Ty"^P(Sl), where 
C^(il) is the set of all infinitely differentiable functions with compact supports in il. We 
use C°^^(ri) and Wi^f{^) to denote the local spaces of functions belonging respectively to 
C°'^(f7') and M^™'P(fJ') for all O' CC Vl. We write u € LP(r2;C) (or u e M^™'P(f7;C), etc.) 
to emphasize that u is a complex valued function. We recall that for m = 1 and p = 2, the 
spaces W^'^(ri;C) and Wq''^{VL]<C), equipped with the inner product 

{u,v) := / S/u-Vv + uvdx, 
Jn 

are Hilbert spaces. Finally, for p > 1 and q being its conjugate exponent, i.e., l/p+ l/g = 1, 
we use W~^''^{rt;C) and Wo~^'''(rJ; C) to respectively denote the dual spaces to WQ'^{il.;C) 
and W^^P{n;C). 

Our main result in this section is the following. 

Theorem 2.1 Let Q C M."^ be a bounded domain. Let 7 G C°'^{T1) for some A G (0, 1). 
Le< N be the Neumann function of L in Q. For y £ CI, denote 70 — "f{y) and let Nq be the 
Neumann function for Lq ~ V ■ 70V ^ ik in CI. Then we have 

\N{x,y)-Noix,y)\<C\x-yf-^+\ Vx e O, x ^ y, (2.2) 

where C is a constant depending only on d,iy,kQ, X,Cl, and ['y]o,\-n- Also, if < \x — y\ < 
dy/2, where dy ~ dist(y, 951), then we have 

I V,(iV(x, y) - No{x, y))\<C {\x - y\'-^+' + k\x - yf-''+^) , (2.3) 

where the constant C depends on diamfi as well. Moreover, if we assume further that 
7 G C^'^{Q), then for all x £ il satisfying < |.t — y| < dy/2, we have 

\V^{N{x,y) ~ No{x,y))\ < C\x - y\'-^+\ (2.4) 
\\7l{N{x, y) - N^{x, y))\<C {\x - y\-''+^ + k\x - y|2-'^+^) , (2.5) 

where C depends only on ||7||ci.^(f2); d,v,kQ,\,Cl, and diamCl. 

In this section, we first consider the Neumann boundary value problems for the operators 
L and its adjoint L* given by 

L* := V • 7V + ik. (2.6) 

Then we give a definition of a Neumann function. Next, we construct Neumann functions, N 
and N* , of respectively L and L* in Q,. Our construction of N and N* holds for a Lipschitz 
bounded domain CI and a coefficient 7 uniformly continuous on CI. If we further assume 
that CI is of class C^, then we are able to derive L^ estimates for the operators L and L* 
with Neumann boundary conditions on dCl. Finally, based on the following global pointwise 
bound for the Neumann function N: 

\N{x,y)\<C\x~y\^~'^ for all x, y G 17 with x- 7^ y, (2.7) 

where C depends only on d^VjCl, fco, and (a modulus of continuity of 7), we describe the 
local behavior of N such as (|2.2p . Assuming that 7 G C^'^{Cl), for < A < 1, we prove that 
estimates ([^ - (^3)) hold. 
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Estimates of (|2.7p -tvDe were derived for the Dirichlet Green's function of L with fc = 
and 7 G L°°(r2) in [TS1[T1]. Under the further assumption that the principal coefficients 
are uniformly continuous of belong to the class VMO, they were generalized to the vectorial 
case in [5J [71 [T31 [H] and to the periodic case in [5J [T^ . 

2.1 Neumann boundary value problem 

We begin with the weak formulation of the Neumann boundary value problem 

J -Lu^ f + V ■ F in n, 

\ {jVu + F)-n = g on dn, ^ ' ' 

where / e Li„^(fi;C), F e ^^^^(l]; C*), and g € Ll^{dn;C). We say that u e Wl^^{n) is a 
weak solution of problem p.Sp if the following identity holds: 



{-l'^u-\/(j) + iku(j))dx^ j {f(j)- F ■'^(j))dx+ j g(j)d(j, V(/) e C°°(f7;C). 
Jn Jon 

Let H ~ VK^'^(17;C). We define the sesquilinear form B{-,-) : H x H ^ C, associated to 
the operator L, as 

B{u,v) := / {'j\7u ■ Vv + ikuv) dx. 
Jn 

It is easy to check that B is bounded and coercive. 

Let / G L2d/(d+2)(^.q^ p g i2(^.(^d)^ g^j^^ g g i2(5^.q^ rpj^gj^ Sobolev 

embedding and the trace theorem, we find that 

T{v) :^ / ifv - F ■VTi)dx + / .giidcr 
Jn Jan 

is a bounded skew-linear functional on H. Therefore, by the Lax-Milgram lemma, we find 
that there exists a unique u (z H such that 

B{u,v)^F{v), yveH. 

We have thus shown that if / e L2d/(rf+2)^^. p ^ L^{n;C'^), and g G_L'^ (dfl; C) , then 
problem (|2.8p has a unique weak solution u in iy^-2(il;C). Since C°°(r2;C) is dense in 
W^'^{il;C), we find that w satisfies following identity: 

( [-i\/u-^ + ikuv)dx= ( [fv- F ■^)dx+ ( gvda, Wv e W^'^{n;C). (2.9) 
Jn Jn Jan 



Let L* be given by (|2.6p . By the same reasoning, we find that there exists a unique weak 
solution w in W^''^{il;C) of problem 

J -L*u ^ f + \7 -F in n, 
j (7VM + F) - n^ g on 

provided / e L2d/(d+2)(-^. ^ ^ ^ L^in^C^), and 5 e L2(ar2;C); i.e., 

{-fVu-Vl>~ikuv)dx^ {fv-F-V^)dx+ gvda, Vv e W^-^{n;C). (2.10) 
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2.2 Definition of the Neumann function 

We say that a function N : fl xil ^ CU {00} is a Neumann function of L in if it satisfies 
the foUowing properties: 

i) N{-,y) e Wl^^iQ) and N{-,y) e W'^'^{n\ Br{y)) for all y € 17 and r > 0. 

ii) N{-, y) is a weak solution of 

f -LN{-,y) = Sy in fl, 



\jVN{-,y) - n^O on dn, 



for all y E ft in the sense 



(7(x)V,iV(x, y) • y(l>ix) + ikN{x, y)(j){x)) dx = ^(y), V</) e C°°(r!; C). 

n 

iii) For any / G C^(r2; C), the function u given by 

u{x) := / N{y,x)f{y)dy (2.11) 
Jo 

is the unique solution in W^''^{^) of problem 

{—L*u = f in J7, 
(2.12) 
"fWu ■ n = on dil. 

We remark that part iii) of the above definition gives the uniqueness of a Neumann 
function. Indeed, let N{x, y) be another function satisfying the above properties. Then by 
the uniqueness of a solution in W^'^{Vl] C) of problem p.l2p . we have 



{N - N){y,x)f{y) dy = 0, V/ e C^{n; C), 

'n 

and thus we conclude that N = N a.e. in Q, x 
2.3 Local boundedness estimates 

Let Bji = Bii{xo) be the ball of radius R centered at Xq, and let u € W^'^(i3fl) be a weak 
solution of —Lu = in Bn. For < p < i?, let 77 be a smooth cut-oflF function satisfying 

< ?7 < 1, suppyycSi?, 77 = 1 on Bp, and \Vt]\ < 2/{R- p). 

By taking rj'^u as a test function, we get 

/ ^rj^\^u\'^ dx ^ - j u ■ V T] dx - ik j Tf\u\^dx. 

J Br J Br J Br 

By taking real parts in the above and using Cauchy's inequality, we get 

/ 777^|Vwp -5R / 2jr]u\/u-\/T]dx <- f 777^! Vwp da; + 2 / j\W r]\'^ \u\'^ dx . 

(2.13) 
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Vu\^dx<-—^^—- I \u\^dx, (2.14) 



Therefore, we obtain Caecioppoli's inequality 

where C — 0(1/). 

Next, we consider the operator Lq defined by 

Lqu ~ V • (70 Vm) — iku = 7oAit — ifcu, (2-15) 

where 70 is a constant satisfying the condition ()2.ip . Let u G W^-'^(B\) be a weak solution 
of —Lqu = 0. Since Lq has constant coefficients, we may apply p.l4p to derivatives of u 
iteratively to get 

< C{m,iy)\\u\\L2(Bi), m=: 1,2,.... 
By the Sobolev embedding theorem, we then have 

sup |m| < C(d)||u||l1/,r.,2(B,/2) < C{d,I^)\\u\\L2(Bi), 

where m = [d/2] + l. Here and throughout this paper [s] denotes the smallest integer not less 
than s. Since the above estimate does not depend on fc, by a scaling argument we conclude 
that if M G W^'^{Bfi) is a weak solution of —Lqu = in Bj^, then we have 

sup |W| <C(d,i/)i?-'^/2||u|U2(B^). 



B 



R/2 



Similarly, if u G W^'^{Bji) is a weak solution of —Lqu = in Bji, then we have 

sup \Vu\ < CR-''/^Vu\\L2(^Bn)- (2-16) 

Br/2 

It follows from the above estimate that for all < p < r < i?, we have 

/ \Vu\'^dx<C{p/rf ( \Vu\^dx, (2.17) 

JBp Jb^ 

where C = C{d,v). Indeed, in the case when p < r/2, we utilize ()2.16p to get the above 
estimate; otherwise, then we may simply take C = 2"^ in (|2.17|) . 

Observe that the same estimates are valid for u G W^'^{B]^) satisfying —LqU = weakly 
in B]^, where Lq is defined as Lq = V • 7qV + ik. 

Lemma 2.2 Assume that 7 S C^{fl) and let 9 be a modulus of continuity of Let Bn = 
Br{xo) C r2 and let u € W^''^{Bfj;<C) be a weak solution of either —Lu ^ f + 'S/ ■ F or 
-L*u ^ f + V ■ F in Br, where f G L1{Br;C) with q > d/2 and F G LViBR/^C^) with 
p > d. Then u is locally Holder continuous in Br and the following estimate holds: 

i?^"Mo,Ao;B„,, < Cq ^^2-d/,||^||^^^^^^ ^ (2.18) 

where Aq G (0,1) and Cq are constants depending on d,iy,p,q, and 9, and [w]o.Ao:n denotes 
the Xo-Holder seminorm of u in D. Moreover, for any po > and < p < R, we have 

sup\u\<c(^iR-p)-''/P» {^j^ |^rdxy^'Vi?^-'^/'?|l/|U,(B«)+i?'-'^/^ll^^llLp(B«)j , 

(2.19) 

where C depends on d,v,p,q,pQ, and 9. 
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Proof. Wc consider the case when u is a weak solution of 

- Lu = f + V -F in Br. (2.20) 

The proof for the other case is identical. Let i?o > be a number to be fixed later. Let 
y G Bi^ and < r < Rq he arbitrary but fixed. Denote 70 = liv) and let Lq be defined as 
in (|2.f 5p . Observe that u is a weak solution of 

-Low-/ + V-F + V-((7-7o)Vu) in Br. 

Let w G Wq {Br{y)) be the unique weak solution of 

f -Lou- = / + V- F + V-((7- 7o)Vu) in Br{y), 
1 w = on dBr{y). 

Then w satisfies the following identity: 

(7o|Vu;p + ifc|-u;p)dx- = / (/W - F • Vw; - (7 - 7o)Vw • Vuj) dx. 

Taking the real parts in the above and using Sobolev embedding, Poincare inequality, and 
Holder's inequalities, wc may deduce that 

7o|Vu|2dx < Cr''^^+'-''^l.f\\LHBAy))\\^M\LHBAy)) 

Br(y) 

+ Cr'''^-'"^F\\L.^BAy))\\^M\LHBAy)) + ^(OI|Vu|U.(b.(,))I|Vu;|U.(b.(,)). 
Denote Ai = 2 — d/q and A2 = 1 — d/p. From the above inequality, wc obtain 

W'^MWABAy)) ^ Cr^'^^^^'\\f\\lABAy))+C^'~^^^^^^\P\\^^^^^ 

On the other hand, observe that v := u — w satisfies —Lqv = weakly in Br{y). Therefore, 
by (|2.f7p . for < /9 < r, we get 



/ \Vufdx<2 \\/v\^dx + 2 \Vw\^dx 

JBAy) JBAy) JBAy) 

<C(-Y [ \Vv\^dx + 2( iVwpdx 
^rl Je^iy) JsAy) 

<c{-\ I \\7ufdx + C [ \\7wfdx 
^rJ JB^(y) JBAy) 

+ ^^"^i X '^•^ + Cr'-^+''^ WfWlABAy)) + Cr'-'^''^~ \\F\ 



LPiBAy))- 

BAy) 

By Campanato's iteration argument (see, for instance, |101 Lemma 2.1, p. 86]), we find that 
if 0{Ro) is small enough, then for all < p < r < i?o we have 

/ \Wu\^dx <C (^) / \Wu\^dx 

JBAy) JBAy) 

I ^ rf-2+2Ao„2(Ai-Ao)|| j-||2 , ^ (i-2+2Ao„2(A2-Ao)|| p||2 

+ r \\J\\LABr{y)) + '-'P ^ 11-^ llLP(B,(y))7 
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where < Aq < min(Ai, A2) = mm{2 — d/q, 1 — d/p). The above estimate (via Morrey's char- 
acterization of Holder continuous functions in terms of Dirichlet integrals; see, for instance, 
POI Theorem 3.5.2]) implies that u is locally Holder continuous in Bji and, in particular, we 
have the estimate 



(2.21) 



Let 77 be a smooth cut-off function satisfying 

< ?7 < 1, suppr/cSfl, r/ = 1 on Bji/2, and iVr]] < 4/i?. 
By taking -q^u as a test function in (I2.20p . we get 



'Br JB__ 



2jriu\/u-\/ridx+ / rj fudx+ ij F-\/udx+ / 27]uF ■ \/ r] dx . 

Br ^ Br ^ Br ^Br 

By taking the real parts in the above and using Cauchy's inequality, we get 

f \\7u\^dx<CR-^ f \u\^dx + CR^ f Ifl^dx + cf \F\^ dx. 

Br/2 ^ Br ^ Br ^ Br 

By Holder's inequality, we then obtain 



I |Vu|2 dx < CR'^ I dx + Ci?2+'i-2<i/9||J||2^^^^^ ^ Ci?'^-2d/p|j^| 

Br/2 ^ Br 



2 

Lp{Br)- 



By combining p.2ip and the above inequality, we get p.lSp via a standard covering argu- 
ment. 

Observe that for any x G 5^/2, we have 

\uix)\ < \uix')\ + \uix)-u{x')\ < |m(x')| +i?^"No,Ao;S«/., V.t' G Br/^. 

By taking average with respect to x' in Bjf/2 and then using (12.181) and Holder's inequality 
we get 



sup 

Br/2 



\u\ < C [r-'^'WuU.^Sr) + R'-^'/'^Wfh.iBR) + R'-'^^FU.^Sr)) 



By using a standard iteration argument (see [11, pp. 80-82]), we obtain (|2.19p from the 
above inequality. This completes the proof. □ 

2.4 Construction of Neumann functions 

The aim of this subsection is to construct Neumann functions of L and L* in and derive 
their basic properties. The following theorem holds. 
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Theorem 2.3 Assume j G C'^^fl). Then there exist Neumann functions N{x,y) and N*{x,y) 
of respectively L and L* inil. Moreover, there exists Xq G (0,1) such thatN(-,y), N*(-,y) G 
^'loc" \ {y}) 0.11 y & ^ and the identity, 



holds. Furthermore, the following estimates hold uniformly in y G fJ, where we denote 
dy = dist(?;, dil): 



i) ll^(-,y)l!L2<</«i-2)(n\B.(j,)) + |iV7V(-,y)|U2(o\B.(^)) <Cri-'^/2 for all r e {0,dy). 

ii) |liV(.,y)|Up(B,_(^)) < Cr^-d+d/P for all r e {0,dy), where p e [I, 

iii) \{x G n : \Nix,y)\ > t}\ < Ct-''/('^-2) fo^ ^ > _ 

iv) ||ViV(-,2/)||iP(B^(.y)) <Cr^-d+d/p for allre{0,dy), where p e [1 , ^) . 

v) \{x<En: \Va:N{x, y)\>t}\< Ct"^'^'^-^^ for all t > dl""^ . 

vi) \N{x,y)\ < C\x - y^^^"^ whenever < \x - y\ < dy/2. 

vii) \N{x,y) - N{x',y)\ < C\x - x'\^°\x ~ y\^-'^-^° if2\x~x'\ < \x ~ y\ <dy/2. 



In the above, C is a constant depending on d,v,ko,i}, and 9; it depends on p as well in ii) 
and iv). The estimates i) - vii) are also valid for N*{x,y). Finally, if q > d/2 and p > d. 



then for any f G L'^{^, C), F G LP{il; C^) and g G Lp'{dVl] <C), the function u given by 



is the unique solution in W^'^{fl) of problem (|2.8p . 

Proof. Wc follow the strategy used in [6], which in turn is based on [13]. Let us fix a 
function $ G C^{R'^) such that $ is supported in Bi(0), < $ < 2, and J^^ dx ^ 1. Let 
?; G be fixed but arbitrary. For any e > 0, we define 



N*{x,y) -.^ N{y,x), 'ix,y n, x ^ y, 



(2.22) 



u 



{x) := / (iV(.T, y)f{y) - VyN{x, y) ■ F{y)) dy + f N{x, y)g{y) da{y) (2.23) 




$,(x) = e-'*$((a;-2/)/e). 



Let v^^y be the unique weak solution in W^ '^{Q,\ C) of problem 



-Lv = <i>e in il, 
7Vw ■ n — on dVl. 



(2.24) 



We define the "averaged Neumann function" N'^{-,y) by 



N\-,y) 



Then N^-^y) satisfies the following identity (c.f. ^^ ): 



I 

Jn 



{j\7N'{-,y) ■ + ikm{-,y)(j)) dx 



I 



nnB,{y) 



^,(j)dx, eW^'^{n;C). 



(2.25) 
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By taking (f) = N''{-,y) = v in (H^Sl), we get 

n Jn JnnB,(y} 

where the last inequahty follows from the Sobolev embedding, namely, 



dx 



Similarly, we get 

k\v\'^dx = 3 /" (7|Vvp + ik\v\'^) dx^'i j ^,vdx < Ce^'^-'^'^/'^\\v\\w^a(n). 

Therefore, we have 

l|iV'(•,2/)llw-.2(0)<Ce(2-'^/^ (2.26) 

where C ~ C{d,u,kQ). 

Let R e (0, dy) be arbitrary, but fixed. Assume that / G C^(r2;C) is supported in 
Br ~ Bfi[y) C ri. Let u be a unique weak solution in W^''^{^l]<C) of problem (|2.12p . We 
then have the following identity (c.f. (|2.10p ): 



{-f'^w ■Vu + ikwu)dx ^ / wfdx, Vw e Ty''^(17;C). (2.27) 
n Jn 

Then by setting = u in (|2.25l) and setting w — iV'(-, y) — v in (|2.27p . we get 



N''{x,y)f{x)dx ^ / (2.28) 
Also, by taking w = m in p.27p . we see that 

/ "f\V u]'^ dx + ik / / ufdx. 

Jn Jn Jn 

Taking the real and imaginary parts in the above and using the Sobolev embedding and 
Holder's inequality 



7|Vu| dx^'St u/dx < C||/||i2d/(d+2)(f2)||u|jwi,2(n), 



k \u\ = 3 / u/rfa; < C||/||i2d/(d+2)(f2)||tt||n/i.2(n)- 
Jn Jn 

Therefore, we obtain 

< C'll/llL2d/(d+2)(f2), (2.29) 
where C ~ C(d, ko). From p.l9|) in Lemma [2^ with po = 2d/ {d — 2), it follows that 

||"||L~(i3«,,) < C {R'-'^'\\uh2.H^-.^n)+R'\\f\\L^(B,)) ■ 
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Furthermore, (j2.29p yields 



||w||i2d/(d-2)(f2) < C R^^'^^^W f\\ I^Br) , 



provided that / is supported iu _B/j. Therefore, by combining the above two inequahties, we 
have 

||«||L~(i3,/,) <Ci?2||/|j^^(^^), (2.30) 

where C depends on d, i/, 17, and 6. By (f2:28l) and (fS^Oj) . we find that for all e € (0,i?/2) 
and R G (0, dy), 



N\-,y)fdx 



Br 



<Ci?2 II /|U»(B«), V/eCr(i3i?;( 



Therefore, by duality, we conclude that 

W{-,y)\\LHBR(y))<CR\ Vee(0,i?/2), Vi?e(0,d,). 

Now, for any x S fi such that Q < \x ~ y\ < dy/2, let us take R 2\x — 2/|/3. Notice 
that if e < R/2, then iV'(-,?/) G W^^^2(5^(^)) aj^^j satisfies ~LN'{-,y) = in Bi?,(a;). Then 
by (|2.19p in Lemma [2.21 we have 

\N^{x,y)\ < Cr-'^||iV^(.,y)|Ui(B.(,)) < Cr-''\\N^{■,y)U^^B,,^^y)) < Cr^-". 

We have thus shown that for any x,y £ il satisfying < |a; — ?;| < dy/2, we have 

\N'{x,y)\<C\x-y\'-^, \fe < \x ~ y\/3. (2.31) 

Next, fix r e {Q,dy/2) and e € (0,r/6). Let 77 be a smooth function on M.'^ satisfying 

< 77 < 1, 77 = 1 on M'' \ Br{y), 77 = on B^/2{y), and |V?7| < 4/r. (2.32) 

We set = = rf-N^{-, y) in (I^T^ to get 

/ 77;^|Vvp dec + ifc / rj^\v\^ dx = ~ I 2^rfuVv ■ Vrj dx. 
Jn Jn Jn 

By taking the real part in the above and using Cauchy's inequality, we get (c.f. ()2.13p ) 

jf]^\\7N%x,y)\^dx <4: I -f\yri\^\N''{x,y)\^ dx. 
n Jn 



We then use (I2.3ip to obtain 

r]^\VN^x,y)\^ dx < Cr-^ 



in J Br{y)\B^/2{y) 

Therefore, for all < e < r/6, we have 

||ViV^(-,y)||i.(o\s.(,))<Cr(2-'^)/2. 



|x-yp(2-'') dxKCr^-". 
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In the case when e > r/6, we obtain from (|2.26|) that 

l|ViV^(-,y)k.(o\s,(,)) < ||V7V^(.,y)|U.(o) < Cr^^-^)/^. 
By combining the above two inequahties, we obtain 

\\^N'{-,y)\\L-in\BAy})<Cr^'''^/'^ Vre(0,d,/2), Ve > 0. (2.33) 
Observe that (|2.3ip also imphes 

\\N^i;y)\\L^.n^-Vin\BAy)) < Cr^''"^/', Ve € (0,r/6). 
On the other hand, if e > r/6, then (j2.26p impUes 

||iV^(.,y)|L../(.-.,(o\s.(,)) < C\\N^i;y)\\w^.^in) < Cr^'-'^y^ . 
By combining the above two estimates, we obtain 

||iV^(-,y)|L=./(.-.,(o\B.(,))<Cr(2'''/', Vre(0,d,/2), Ve > 0. (2.34) 

From the obvious fact that dy/2 and dy are comparable to each other, we find by p.33p and 
(|2.34p that for all < r < dy and e > 0, we have 

||iV^(-,y)|L../(.-.,(a\B.(,)) + \\^N^i-,y)\\LmABAy)) < CA^~''^/\ (2.35) 
From (p:^ it follows that (see [H PP- 147-148]) 

\{xen:\N'ix,y)\>t}\<Ct-^^^'^-^\ Vi > Ve > 0, (2.36) 

\{xen:\S/^N'{x,y)\>t}\<Ct-'^^^'^-^\ Vi > Ve > 0. (2.37) 

It is routine to derive the following strong type estimates from the above weak type estimates 
([Oell and ([237)) (see, for instance, [El p. 148]): 

\\N'{;y)\\L.iBAy)} < Ct'-'+'^p, Vr G (0,rf,), Ve > 0, Vp £ [1,^), (2.38) 
l|VA^'(-,2/)||iP(B.te)) <Cri-'^+'^/P, Vre (0,dy), Ve > 0, Vp G [1, /y). (2.39) 

From (|2.33l) . (j2.38p . and (|2.39p . it follows that there exists a sequence {e„}^]^ tending to 
zero and a function N{-,y) such that N'^"{-,y) converges to N{-,y) weakly in W^''^{Br{y)) 
for 1 < p < d/{d — 1) and all r G (0, dy) and also that A^*^" (•, y) converges to A^(-, y) weakly 
in iy^'^(r2 \ Br{y)) for aU r e (0,(iy); see [13l p. 159] for the details. Then it is routine to 
check that N{-,y) satisfies the properties i) and ii) at the beginning of Section [2?2| and also 
the estimates i) - v) in the theorem; see [131 Section 4.1]. 

We now turn to the pointwise bound for N(x, y). For any x Cz Q such that < \x — y\ < 
dy/2, set R 2|a:;-y|/3. Notice that ((235)) implies that N{-,y) e W^'^{Bh{x)) and satisfies 
—LN{-,y) = weakly in Bfj{x). Then, by p.l9p in Lemma [2.21 and the estimate ii) in the 
theorem, we have 

\N{x,y)\ < CR-''\\N{;y)h.^Bn{cc)) < CR-'\\N{;y)U.^Bs,iy)) < C|.t - yp- f 

We have thus shown that the estimate vi) in the theorem holds. Then, it is routine to 
see that the estimate vii) in the theorem follows from p.lSp in Lemma [2.21 and the above 
estimate. 
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Next, let X e r2 \ {y} be fixed but arbitrary, and let N" (-,2;) € W^''^{n;C) be the 
averaged Neumann function of the adjoint operator L* in f2, where < e' < d^. Then we 
have 



{-fVN' iz,x)-Vi}{z)-ikN'' {z,x)'iP{z))dz= $,'(z)V;(z)d2, (2.40) 

JnnB^, (an) 

for ah ip G W^i'2(fJ;C). By setting 7V^'(-,x) in (P?^ and -0 = N''{-,y) in (piigi) and 
then taking complex conjugate, we obtain 

[ ^^:N''{-,y)dz ^ [ '^>^m'{-,x)dz. 

Let N*{-,x) be a Neumann function of L* in f2 obtained from A^'^™(-,x), where {em.}m=i is 
a sequence tending to 0. Then, by following the same steps as in [131 P- 151], we conclude 



N{x,y)^ N*{y,x), 

which obviously implies the identity (|2.22p . We remark that by following similar lines of 
reasoning as in |13[ p. 151], we find 

N^ix, y) = ^ (^) Nix, z) dz, 

and thus we have in fact the following pointwise convergence: 

lira N^{x,y) ^ N{x,y), Vx, y e 17, x^y. (2.41) 

Now, let u be the unique solution in W^^'^{Vl;C) of problem (^1^ with / e C;?°(f7;C). 
By Lemma [2.2[ we find that u is continuous in fi. By setting w = N'^{-,y) in (|2.27p and 
setting (/) = 7i in (|2.25p . we get 



N^{x,y)f{x)dx— j ^eudx. 
We take the limit e — > above and then take complex conjugate to get 

u{y) = / N{x,y)f{x)dx, 



which is equivalent to p. lip . We have shown that N{x,y) satisfies the property iii) in 
Section [2.2| and thus that N{x, y) is the unique Neumann function of the operator L in fl. 

Finally, let / S L''{Q,] C) with q > d/2 and g S L^{di}; C), and let u be the unique weak 
solution in W'^-'^{il;C) of problem ((Zl)) : see Section[2TT] Then u satisfies the identity ([2^ . 
By setting v = N'^ (•, a^) in (|2.9p and setting ijj ^ u in (|2.40p . we get 

(N^'{z,x)~WN'''{z,x))f{z)dz+ N^'{z,x)g{z)d(T{z)^ 1 ^^.udz. 
^ ^ Jdn JnnB^,{x} 

By Lemma l2.2[ we again find that u is Holder continuous in fl. Then by proceeding similarly 
as above and using (|2.22p . we obtain 



u{x)^ / iN{x,y)fiy) - VyNix^y) ■ F{y)) dy + / Nix,y)giy) da^y), 
Jn Jn 

which is the formula (|2.23p . The proof is complete. □ 
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2.5 LP estimates 

We now assume that is a bounded domain. In the foUowing lemma we obtain 
estimates for the operator L with uniformly continuous coefficient 7. 

Lemma 2.4 Letfl C Mf^ be a houndedC^ domain and assume that'y S C^{rt). Let q G 

p G (1,00), and s = min(g*,p), where q* = dq/{d— q). For each f G L'^(f2;C) and 

F £ LP(f2;C'*), there is a unique weak solution u G W^''^{Vl) to 



-Lu = f + W ■ F in n, 
(7VU + F)-n = on dn. 



(2.42) 



Moreover, the foUowing estimate holds: 

lkllwi.»(o)<C(||/|U,(o) + ||FllLp(o)), (2.43) 
where C depends on d, v, ko,p, q, fi, and 9. 

Proof. Note that in the case when / = 0, the proof for estimate (|2.43p reduces to 

hlk-i.p(o) <C\\F\\L.in)- (2.44) 

In this case the proof for the existence and uniqueness of weak solution u G W^''P{fl) as well 
as the estimate (|2.44p follow essentially from the same argument as in [17] . 

We consider the case when / is not identically zero. Observe that L'?(r2) C i^) 
with the estimate 

||/||vv-i„*(a) <C||/||L.(n), where C^C'id.n). 

Then by [51 Corollary 9.3], there exists a unique weak solution v in W^'"^ (il) of the Neumann 
problem 

1 



dv/dn = on d^l, 
where \^\ is the volume of fl. Moreover, v satisfies the estimate 

l|VHL,.(o) <C||/- <C||/|U,(o). (2.45) 

Then, we apply estimate (|2.44p with F + Vw + / dy)x and s in place of F and p, 

respectively, and use Holder's inequality to get estimate ()2.43p . □ 

We denote by LP'°°{fl) the usual weak space. The following lemma is a variant of 
Lemma \TM in the weak Lebesgue spaces. 

Lemma 2.5 Let C be a bounded domain and assume that 7 G C'^{n). Let q G 
{l,d), p G (1,00), and s = mm{q*,p), where q* ~ dq/{d — q). If f £ L'''°°{^l;C) and 
F G LP'°°(fi; C*), there is a weak solution u of problem (|2.42p that satisfies an estimate 

||Vu||l.,^(o) < C (||/||L9-=e(o) + ||F||ip,oo(f2)) , 
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and, for s < d, the following estimate as well: 

\\u\\L''*-'>°{n) < C" (i|/||L9.°=(o) + \\F\\Lp-°^{n)) ■ 

Moreover, there is uniqueness of weak solutions to (j2.42p in the sense that if u is a solution 
in W^'*(fl) for some t > 1, then u = u. 

Proof. The lemma follows immediately from Lemma 12.41 by applying [7j Lemma 1] to the 
solution operator T : F u as well as to the map / h-> u in (|2.45p . □ 

Lemma 2.6 Let U, and 7 satisfy the same assumptions as in Lemma \2.4\ There exists a 
constant Ci > such that the following holds: For any f G C^(fl; C), let u € W^''^{fl; C) be 
the unique weak solution of 

^Lu ~ f in Q. 1 —L*u = f in fl 

or < 

7VU • 71 = on ofl I 7VU ■ n = on ofl. 
Then for all x ft and < R < diam(ri) , we have 




The constant C'l depends on d, ly, fl, and 9. 

Proof. We will only consider the case when w is a weak solution of —Lu = f with zero 
conormal data. By Lemma 12.41 we find that u G W'^'P{n) for ah p e (l,oo) and 

l|Vii||iP(s-2) < C||/||ipd/(p+d)(f^) < C||/||ioo(Q). 

Let V = (u, where ^ : R"^ — > M is a smooth function to be chosen later. Observe that is a 
weak solution of the problem 

-Lv = / + V • F in n, 
(7VW + F) • n = on dn, 

where 

/ := Cf - 7VC • Vm, F := -7^/VC. 

Let X G and < R < diam(O) be arbitrary but fixed. For any y E O Bnix) and 
< p < r < R, wc choose the function C to be such that 

0<C<1, suppCcB,(y), C^l on Bp{y), and |VC|<2/(r-p). 

For any p E {l,oo), we set q — pd/{p + d) and apply Lemma [2.41 together with Holder's 
inequality to get 

||Vu|Up(n,) < C (ri+'^/ni/||L~(o.) + {r - p)"'!! V7i|Up./(p+.,(n.) + (r - p)-'\\u\\L.in)) , 

(2.46) 
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where we use the notation fir ~ ^riv) = H Br{y). Now, Gx p > d, and let m 
[d(l/2-l/p)], 

Pd , , {r - p)j . 

= — : and rj=pH , j^O,...,m. 

a + pj m 

Then we apply (|2.46p iteratively to get 

Notice that 1 < p,„ < 2. By using Holder's inequality we then obtain 

/ \ m — 1 

If we take r = R/A and p < r/2 = R/i in the above, then for all y e 7^/4(0;), we get 



s:2p(y) 

+ + Ci?^(i/^-i/2)||Vii|U.(o,,,(.)) =: (2.47) 

Hereafter in the proof, we shall denote fi/j = Ofl(x). Then by Morrey-Campanato's theorem 
(see [m Section 3.1]), for all z, z' G we have 

\uiz)-uiz')\ <CR^-'^'PA{R), 

where A{R) is as defined in (|2.47p . Therefore, for any z G 51^/4 wc have 

\u{z)\ < \u{z')\ + \u{z) - u{z')\ < \uiz')\ + CR^-'^'PA{R), Vz' e ^r,^. 

By taking average over z' e ^r/a in the above and using the definition of A{R), we obtain 

sup 1^.1 < / |^.(z')l^^2' + Ci?'ll./l|L-(f^.)+Ci?-''/nk||LP(I2,)+Ci?l-'^/2||V7.|U.(n,,,). 

Then by using Holder's inequality and Caccioppoli's inequality, we get 

sup \u\ < CR'WfWL^^nn) + CR-''^P\\u\\L.(nn,) + CR-''/^ 

By using a standard argument (see |1H pp. 80-82]), we derive from the above inequality 

sup \u\<CR^\\f\\L^^nR) + CR-'/-'\\ 

The proof is complete. □ 



16 



11 
iv 



2.6 Global estimates for Neumann function 

The next theorem provides global pointwise bound for the Neumann function N. 

Theorem 2.7 Let il CZ M."^ be a bounded domain and assume that 7 G C^{il). Let 
N{x,y) be the Neumann function of L in Q as constructed in Theorem \2.3[ Then we have 
the following global pointwise bound for the Neumann function: 

\N{x,y)\<C\x-y\^-'^ for all x,y e n with x ^ y, (2.48) 

where C depends on d, v, fi, and 9. Moreover, for all y ^ and < r < diam(r2), we have 

i) ||iV(.,y)|L2./(.-2,(o\s.(,)) + \\^N{-,y)\\L2(n\BAy)) < Cr^-"'^ . 

ii) ||A^(-,2/)|Up(onB.(,)) <Cr2-^+''/f forpe[l,^). 
i) \{xen: \N{x, y)\ > t}\ < Ct-'^/t''-^) for all t > 0. 
/) ||ViV(-,2/)|Up(nns.(^)) < Cr^-d+d/P forpe[l,^). 

v) \{x e n : \V^N{x,y)\ > t}\ < Ct-rf/('i-i) for all t > 0. 

vi) \N{x,y) - N{x' ,y)\ < C\x - x'\^°\x - yp-^-Ao j^, _ ^'^ < \x - y\/2 for some 
Ao e (0,1). 

In the above, C is a constant depending on d, v, fco, fi, and 9; it depends on p as well in ii) 
and iv). Estimates i) - vi) are also valid for the Neumann function N*(x,y) of the adjoint 
L*. 

Proof. Let y S be arbitrary, but fixed. Assume that / G C^{fl;C) is supported in 
^niy) = 5^ n Bji{y) and let u be the unique weak solution in W^ '^{Q.] C) of problem (|2.12p . 
Then we have the identities (|2.27p and (|2.28p as in the proof of Theorem l2.3l Also, we have 
estimate (|2.29p . and thus by Sobolev embedding theorem, we get 

< C!|/|L../<.+.,(f,) < (2.49) 

where C = C{d, v, fi). Then by Lemma [2.61 and (|2.49p . we obtain 

l|w||L-(f2„/2(y)) <CR^\\f\\L--(na(y))- (2.50) 

Hence, by ((^^ and we conclude that 



/ N'[x,y)f{x)dx 



< CR^WfU^^n^Xy)), V/ e C^{nn{y);C), Ve G (0,i?/2). 

(2.51) 



Therefore, by duality, we conclude from (j2.5ip that 

\\N^{;y)\\LHnniy)}<CR^, Vee(0,i?/2). (2.52) 

Next, recall that the v = N'^{-,y) is the unique weak solution in W^''^{rt;'C) of problem 
([2211). Let X e n, r > 0, and e > be such that B^{y) n Br{x) = 0. Then Lemma \T6\ 
implies that 

l|A^^(-, y)l|L~(o.,,(.)) < Cr-''^'\\N^{; y)\\LHnAx)y (2.53) 
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By a standard iteration argument (see [11] pp. 80-82]), we then obtain from (|2.53p that 

\\N'{-,y)\\L^(n^M-)) ^ C'^"'ll^'(-.y)IU^(fV(x))- (2-54) 

Now, for any X en \ {y}, take R = 3r = 3\x - y\/2. Then by ([^31)1 and ^Ml, we obtain 

for aU e G (0, r) that 

|A^^(a;,y)| < Cr-'*||iV'(.,y)|Ui(o,(,)) < Cr-'^||iV^(.,2;)|Ui(n3.(,)) < 

Therefore, by using (|2.4ip . we may take the hmit e in the above and obtain (|2.48p . 

To derive estimates i) - vi) in the theorem, we need to repeat some steps in the proof 
of Theorem 12. 3l with a httle modification. Let v ~ N'^{-,y), where < e < mm{dy,r)/6 and 
< r < diam(r2). Let be a smooth function on R"^ satisfying the conditions (|2.32p . We 
set 4) = rfv in (|2.25p and obtain 

/ (777^ Vw • Vw + ikvv) dx + 2ri^vVv ■ Vrj dx = 0, 
Jn Jn 

where we used the fact that rj^^^ = 0. By using Cauchy's inequahty we get 

/ ij'\yN'{;y)\'dx<C f \V,f\N'{;y)\Ux. 
Jn Jn 

By using the pointwisc bound for N'^(x,y) obtained above, we get 

/ \VN%-,y)\^dx<Cr-^ f \x - yl^-^'^ dx < Cr^-'^ . 

Jn\BAy) JBr{v)\B^/2{v) 

By taking the hmit e — > in the above, we get 

\\'^N{-,y)\\mn\BAv)) < Cr^^-'^)/^, < Vr < diam(f7). 

Observe that the pointwisc bound (|2.48p together with the above estimate yields 

\\N{■,y)\\L2.,i.-2^n\BAv)) + Ni;y)\\L2^n\BAv}} < Cr'-^^''^/^ < Vr < diam(17), 

(2.55) 

where C depends on d, v, f2, and 9. 

By following literally the same steps used in deriving p.36p - (|2.39p from (|2.35p . and 
using the fact that \n\ < 00, we obtain estimates i) - v) from (|2.48p and (|2.55p . 

Finally, we remark that the proof of Lemma l2.6l in fact implies that there exist constants 
Aq G (0, 1] and Ci > 0, which depend on d, fi, and 6', such that for all x £ 51 and 
< i? < diam(il), the following holds: Let u be a weak solution in VK^^^(ilfl'(x)) of cither 

-Lu = in finSfl(x), 7Vu-n = on dfinBii{x) 
or -L*u = in CinBu^x), 7VU • n = on dCinBii{x), 

then wc have 

By utilizing the above estimate and modifying the proof for estimate vii) in Theorem 12. 3| 
we have vi), and the proof is complete. □ 
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2.7 Proof of Theorem 

We are now ready to prove Theorem 12. II 

Let u = N{-,y) - No{-,y). Observe that Theorem [Q imphes that u e W^^'^ifl) for 
1 < q < d/{d — 1), and also that we have 

/ {jS/uV$ + iku$)dx = I (7-7o)V7Vo(-,y)Wrf2:, V0€ C°°(f];C). 
Jn Jn 

In other words, m is a weak solution in W^''^{Vl) of the problem 

-Lu = -V ■ F in n, 
(7VU + F) • n = on on, 

where F = (7 - 7o)ViVo(-, y)- 
Note that 

\y.No{x,y)\<C\x~y\'-'', Wx^yeQ, x^y. (2.56) 

Indeed, for any a; <E SI with x ^ y, we set R= \x — y\/2 and apply (|2.16p and estimate i) in 
Theorem 12.71 to obtain 

|V,7Vo(a;,y)| < CR-^/'^WV No{-,y)\\mn\B,,(y)) < CR'-", 

which obviously implies (|2.56p . Moreover, by repeating the same argument, we have 

\V'^No{x,y)\<C\x-y\^-''-\ Va;,yGl}, x^y, fc = l,2,.... (2.57) 

We then obtain 

\Fix)\<C[-f]o,x-n\x-y\-''/", Vx G fi, x^y, (2.58) 

where a = d/{d— 1 — A), and hence F e L'^{Vl) for all q < a. It then follows from Lemma [2^ 
that u £ W^''^{Vl) for all q S (1, a). In fact, by Lemma [2.51 we have 

hllL=*-(o) + l|Vu|Uo.~(o) <C. (2.59) 

Let V = Cu, where C : I^'' — > R is a smooth function to be fixed later. Observe that v is 
a weak solution of the problem 

-Lv = f + V -F in n, 
(7VU + F) ■n = on dn, 

where 

/ -VC-F-7VC- Vu, F CF --fuVC- 
Notice that if C = on a neighborhood of y, then we have / £ L'(ri) and F G L'' (SI) for all 
q e {l,a). By Lemma we have v e W^'"^' (ft) and thus, we find that u £ W^^l (fi \ {y}). 
By repeating the above argument, if necessary, we conclude that u £ Wl^^{^\{y}) for some 
s > d, and thus we have u £ L1^^{Q \ {y}). 

Next, for x £Q with x ^ y,\ct R^ \x-y\/2. For any x' £ nr]B]i{x) and < p < r < R, 
we choose the function C to be such that 

0<C<1, suppC C B^(a;'), C=l on Bp{x'), and |VC| < 2/(r - p). 
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Notice that for all q e (1, d), we have the following estimates, where we write Vlp = Vlp{x') = 
n Bp{x') for the simplicity of notation, 

I|Vu||l<!*.-(j2j < ||Vu||i,*,o=(n), 
||VC-F|U„^(o) < ||VC||L~||i^|U„~(n.) <C(r-p)-ir^/«|lF|U^(a.), 
live- Vm|U,,»(o) < llVC|lL-|iV?/|U,,.^(fi^) < C(r-p)-i||Vu|U,,.^(n,), 

||uVC||l,-,oo(o) < ||VC||L-||u|lL,-,~(n^) < C(r - p)"i||it||i<,..oo(Q^). 
Therefore, by Lemma [^31 applied to u, we have for all t G (d/ {d — 1), od) 

l|V7.|U*,^(o,) <c{{r- p)-Vi+'^/'||F|U^(j,,_) + (r - Vu|L../(.+.,,^(^,_) 

+r-''/*||^^||L-(o.) + {r- p)-i||u|U*,oo(j,^)) . (2.60) 
Now, fix s > d and let m = — 1/s)], 

= '^^ . and rj = p + j = 0, . . . , m. 

a + sj m 

Recall that if is a bounded set and < g < p < oo, then 
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HE) < W/Wlhe) < y^l^l'^'"'^' ll/IU..~(i^). (2.61) 



With the aid of (|2.6ip . we apply (|2.60p repeatedly and argue as in the proof of Lemma 
to obtain 



/ \ m / \ m — 1 

If we take r = R/A and p < r/2 ~ R/A in the above, then for all x' £ f2j^/4(x), we get 

+ +Ci?'^(i/^-i/")||Vu|U..oo(o), (2.62) 

which is analogous to (|2.47l) in the proof of Lemma 12.61 Then by utilizing (|2.59p and 
proceeding as in the proof of Lemma 12.61 we obtain 



sup \u\ < CR-^'WuWmnu) + CR\\F\\L^^^nu) + CR^-'"'^, ^ ^nix). 
By Lemma [2.51 and (I2.6ip again, we get 
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Combining the above two inequalities and using (|2.58p . we get 

\N{x,y)~No{x,y)\ ^ \uix)\ < CR^^'^+^ + CR^-^''^ < C\x ^ y\^~'^+^ . 



This completes the proof of (|2.2p . 

Next, we turn to the proof of (|2.3p . Let u = N{-, y) — Nq{-, y) as before. Observe that u 
satisfies 

~Lqu V • (F + (7 - 7o)V-u) in Q. 

Let X G n satisfy < \x — y\ < dy/2 and let R ^ \x — y\/2 as before. For any x' £ Bii/2{x) 
and < r < R/2, let w be the imique weak solution in Wq'"^ {Br{x')) of the problem 

J -7oAw = -i/sit + V • (i^ + (7 - 7o)Vu) in Br[x'), 
\ w = on dBr{x). 

Then w satisfies the following identity: 



70 Vw • Vw dz = - / (ikuw - [F - Fr) • Vw - (7 - 7o)Vu • Vw) dz, (2.63) 



iBAx') 

where we use the notation 



F — F 

J. r - J. .j; 



\Brix') 



Fdz. 



Br{x') 



Notice that by Holder's inequality and the Sobolev inequality, we have 



L 



ikuw 



Br{x') 



< Ck 



(d+2)/2d 



^|2d/(d 



+2) 



IB Ax') 

<CkA'+'y^uU^^s,) 
Also, by Holder's inequality, we have 



1/2 



\Vw\ 



BAx') 

\ 1/2 
2 1 



BAx') 



{F ~ Fr) ■ Vw 



BAx') 



< 



1/2 



BAx') 



\F-Fr\ 



1/2 



|Vu;|- 



BAx') 



<[F]o,X:Bs.rW/'( f IV 

\JbAx') 



1/2 



1/2 



BAx') 



\Vw[ 



Similarly, we estimate 

(7 ~ 7o)Vu • Vw 

BAx') 



< II7- 7o||l~(b,{^')) 

< b]o,X;B^r^ I I 

\JbAx') 



1/2 



BAx') 



\Vu\ 



1/2 



BAx') 



\Vw\ 



1/2 



1/2 
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Therefore, by using Cauchy's inequalities, we derive from (j2.63p and the above estimates 
that 

(2.64) 

where we use abbreviation Bb. = Bb.{x). 
Notice that v = u — w satisfies 

Av = in Br{x'). 

By well-known estimates for harmonic functions (see, for instance, |10[ p. 78]), we get 
\Vv-iVv)p\^dz<C{p/rf+^ [ \Vv~ {Vv)rfdz, Vpe(0,r). 

B^{x') Jb^x') 

Then by using the triangle inequality, we get for all < p < r that 

/ |Vu- (Vw)pp(iz < 2 /" \\/v-{\7v)p\^dz + 2 [ \\/w-{\/w)p\'^dz 

JBp(x') JBp(x') -JB^ix') 

<C{p/rf+^ [ \Vv ~ {Vv)r\^ dz + 2 [ \Vw\^dz 

JbJx') JbJx') 



< C'{p/r) 



d+2 



[ \\7u- {\7u)r\^ dz + C [ iVwpdz, 

Jb„(x') J Br{x') 



o{x') JBr{x') 

where we have used the well known fact that 



inf / |/-c|2dz= / \f-fr\^dz. 
"^^JbAx) JbAx) 



r{x) 

By combining the above inequality and ()2.64p . we get for all < p < r that 

\Vu-{Vu)p\'^dz<c{-Y^^ [ \yu-{\/u)r\^dz 
Bp{x') ^rJ JbAx') 

JbAx') 



On the other hand, by setting e = d/s and p = r in p.62p . we get 

1/2 



f \Vufdz\ <CR^\FUo.^B„)+CR'-'\\u\\L^^Bn) + CR'-''^". 

JB,.ix') J 

Combining the above inequalities, for all x' e -8^/2(3;) and < p < r < R/2, we get 

/ \yu-{S7u)p\^ <c(^y^' [ |V^-(V^.),f+CfcV''+2h|li.(5^)+C[i^]L.B. 

■'Bpix') Jb^x') 

+ C[7]L;S„ r'+''-'' + R^-'\\u\\^^^B.^ + R^~'''-' 
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By Campanato's iteration lemma, for all x' G Bjf/2{^) a-nd < ?' < -R/2, we have 
where we set /? := A — e € (0, 1). Therefore, by Campanato's theorem, we obtain 

R^[Wu]o.,(i,B^,, < Ci?-'^/2||Vu|U2(B^) + CkR\\u\\L..iBn) + CR\F]o,X,Bh 

By Caccioppoli's inequality, we estimate 

||Vu|U.(B„,) < CR-'\\uU.^B,,,^_)+C\\FU.^B,„,^ < CR'''^R'-''+\ 
Also, observe that 

[F]o,X;Bn < i?^[7]o,A;B„[V7Vo(-,y)]o,A;B, + Ho, A;B« || ViVo(-, 2/) < Ci?!-'*, 

where C depends on d^v^X^Q, and [7)0, A;0- Therefore. 

i?^[Vu]o,/3;BH/2 < C {R^-'^+^ + A:i?3-'^+^ + < + fci?^), (2.66) 

where we used the assumption that Q is bounded in the last step. By proceeding as in the 
proof of Lemma 12.61 we derive from (|2.66p that 

sup |Vu| < CR^~'^+^{1 + kR^). 

Br/4, 

This completes the proof of (|2.3p . 

Now, let us assume that 7 e C^''^(n). Let x E satisfy < \x ~ y\ < dy/2. We again 
set R = \x — y\/2 and write = Bb{x). Observe that u satisfies 

-7AU = / in _B_R, 

where 

/ := V7 • Vit - iku + ik%^{-f - 7o)iVo(-, y) + V7 • VNo{-,y). 

We claim that / € C^'^{Bp). Indeed, observe that by feeding estimate (|2.3I) back to p.65p 
and repeating the above steps, we obtain an improved version of estimate ()2.66p . namely, 

[VM]o,A;i3«/, < CR^-\l + kR^). 

Therefore, we obtain 

[V7 • Vu]o,A;Sb < [V7]o,A;I2||Vm||loc(b^) + ||V7||L^(n)[Vu]o,A;BH 

< Ci?i-''+^(l + kR^) + CR^-'^{1 + kR^) < CR^-'\l + kR^), 
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where we have used the assumption that Q is bounded. Also, by taking s ~ d/{l — X) in 
(|2.62p . we find that for all x' G -B_r/2 and p < i?/4, we have 

/ -(d-2+2A) f l^^l^A <Ci?l-^||F|Uoo(s,)+Ci?-^h|U=o(B,)+Ci?2-'i<Ci?2-'^. 

From the above inequality and a standard covering argument, we find that 

[iku]o,X;Ba = k[u]o^x;BR < C'kR^~'^. 
In a similar fashion, with the aid of (j2.57p . we also estimate 

[ikio'h-lo)Noi-,y)]oA:B^ < CkR^-", 
[V7-ViVo(-,2/)]o,A;S« <Ci?i-'^-\ 

Combining all together, we find 

[/]o,A;B„ < CR^-'^{1 + kR^) + CA:i?2-d + CR^"^'^ < CR-'^{\ + fci?^), 

where we again used that diamfi < oo. Then the interior Schauder estimate yields 

[V2u]o,A;B,/, < C ([/]o,A;B„ + i?"'"^ h|Uoo (s^)) < CR-\l + kR^). 

On the other hand, by the standard estimates, we have 
W^^Al^b^,,) < C (i?-i||V«|U2(B„) + ||V«|U2(B„) + ||VF|U2(s„)) 

< C + R-')\\u\\mB:,,,,) + + l)!|i^!U=(B3H/.) + ||VF|U2(B„)) 

< CR'^^^R-'^+^il + R^-^ + R)< CR'^/^R-'^+^. 
Therefore, we have 

sup \V^u\ < CR-''/^\\VMl^(b^,,) + CR^[V^ 

Br/4, 

< CR-'^+^ + CR-'^+^{l + kR^) < CR-'^+^il + kR^). 

We have thus proved (|2.5p . Finally, we prove (|2.4p as follows. Notice that v :— du/dxi, for 
i = 1, . . . ,d, satisfies 

-Lv = V • F, where F = {dj/dxi)Wu + dF/dx^. 

Let R= \x — y\/2 as before and applying (j2.19|) in Lemma [2T2l to v, we obtain 

sup \v\ < C (i?-'^/2||Vzi|U2(B,) + R'-'^'\\FU2^Bn)) ■ 
Br/2 ^ ' 

Notice that 

\\F\\L-iBR) < CWV^h^^n) (l|Vii|U2(B^) + \\S/Na{-,y)\\mBR) + R\\^'No{-,y)\\L^^BR)) 
<C||Vu|U2(s,)+Ci?'^/'i?i-'^. 
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On the other hand, by (|2.62p . we find 



Combining together, we obtain 

sup |Vm| < C {R^-'^+^ + i?2-'i+A ^ < Ci?^-''+^, 

Br/2 

where C depends on |j V7||Loo(f2) and diamil as well as on d, v, fcg. A, f2. This proves estimate 
(j2.4p . The proof of Theorem 12. II is now complete. 

Remark 2.8 We remark that for z G O fixed, we may choose e = e{z) > so small that 
for all y G B^{z), we have dy > 4e. Then all x,y G B^{z) should satisfy the relation 
\x — y\ < dy/2, and hence i2.2\) - ![275\) hold for all x, y G Bf^(z). We also note that in the 
proof of ()2.4p . it is enough to assume that 7 G C^(f2) not 7 G C^'^{^1). Also, if we assume 
7 G C'^{fl), then instead of (|2.5|) . we have 

\VlN{x,y)~VlNQ{x,y)\ <C\x-y\-'^+^, Va; G 17 satisfying < \x -y\ < dy/2, (2.67) 

where C depends on \\^\\c'^(n) and diamfi as well as on d,v,ko,\,ri. The proof for (j2.67p 
is analogous to that for \2A\ . Moreover, if d > A, then we may take 7 = 1 in p.2p . p.4p . 
and (|2.67p since in that case, we may take a = d/ {d — 2) < d in (|2.59p . 



3 Applications to quantitative photo-acoustic imaging 

In this section we deal with the problem of quantitative photo-acoustic imaging to recon- 
struct the optical absorption coefficient from the absorbed energy density. The absorbed 
energy density can be reconstructed using the measurements of the acoustic wave on the 
boundary of the medium. See, for instance, [II I21j. 

Reconstruction of the optical absorption coefficient, /x^, from the absorbed energy density, 
A{x), is subtle since fia is related to Aix) in a nonlinear and implicit way. In fact, /^a is 
related to A{x) by 

A = /ia$ (3.1) 

Here $ is the light fiuence which depends on the distribution of scattering and absorption 
within f2, as well as the light sources. Let fis be the scattering coefficient. The function $ 
is related to /ia through the diffusion equation 

T+^°(")-^^- M(x)L(.) ^)''("^=° 
^ o l^a\X ) ^ ^sydj) J 

with the boundary condition 

1 9$ 

5 on d^, (3.3) 



3(/iQ(a;) -I- [is{x)) dv 



where g denotes the light source and w is a given frequency. Equation p.2p is derived 
based on the diffusion approximation to the transport equation which holds when fi^ ^ fJ-a- 
Sec, for instance, [H [M]. Note that in [2], the boundary condition is a Robin boundary 
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condition. However, it is easy to check that ah the estimates derived in [5J Section 2] hold 
for the Neumann boundary condition p.3|) . 

We restrict ourselves to the three-dimensional case and suppose that the medium contains 
a small absorbing anomaly whose absorption coefBcient is to be reconstructed. The small 
unknown anomaly D is modeled as 

D = z + eB, (3.4) 

where z represents the location of D, B is a reference domain which contains the origin, 
and e is a small parameter representing the diameter of the anomaly. We assume that the 
anomaly is away from the boundary dil, namely 

dist(z, an) > Co (3.5) 

for some constant Cq. Since D is small and absorbing, and the background absorption is 
quite small compared to the scattering, we may assume that 

fia{x) = fiaXoix) (3.6) 

where /ia is a constant and xd is the characteristic function of D. Then, p.2p and p.3p 
may be approximated by 

(—+ HaXDix)-lv ■ ^^v]<^{x) =0 inO, 

1 9$ ^-^-'^ 



g on dfl. 



3/Xs(a;) di^ 

Since D is small, we may regard $ as a perturbation of which is the solution of 

!f^^„iv^V^$W(a:) =0 inn, 

" ^'^"^ ^ ,0. (3-8) 

- g on ail. 



3^six) dv 



The reconstruction methods in [2] deeply rely on the following asymptotic formula $ — 
which was obtained under the assumption that /is is constant: 

($ _ ci>(0))(2) ^ 3e2/i„/i,$(o)(z)iVB(0) - e^SsWm ' V$(o)(z), (3.9) 

Ms 

where Nb be the Newtonian potential of _B, which is given by 

Nb{x):= / T{x-y)dy, a; e K^, (3.10) 



and Sb is the single layer potential associated to B, which is given for a density ijj £ L^{dB) 

by 

SBm{x):^[ T{x-y)^y)da{y), x 

JOB 
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Here r is the fundamental solution to the Laplacian in three dimensions, i.e.. 

1 



4:Tt\x\ 



The purpose of this section is to show that the asymptotic expansion (|3.9p holds even 
when fig is variable. The following theorem holds. 

Theorem 3.1 Let be a bounded -domain in R"^. Let D ~ z -\- eB, where B is a 
bounded Lipschitz domain in M'^ containing the origin. Suppose that *s given by iS. 6\) 
and Us S C^''^(f2), A G (0, 1), and set jig l^s{z). Then, as e — > 0, we have 

($ _ $(0))(2) « 3eVaMs*^°H^)^s(0) - e^SeMiO) ■ V^^°\z), (3.11) 

Ms 

where the error term is less than 

aMs + ( — ) + e Ma 
Ms / 

for p > 3 and some constant C depending on ||ms ||ci.^ j -^j ^/c, and g. 

Since the proof is essentially the same as that in [5], we only outline the proof without 
much details. 

Let N{x, y) be the Neumann function of the operator ^ — i V • V on Vt. Then one 
can show by following the same lines of the proof of [H Lemma 2.1] that for any cc € 51, 

($-$W)(x) = Ma / ^{y)N{x,y)dy 

If 1 1 (3-12) 



(2 - , Ma\ 




£ ^iat^s + — 


+ e^/JTsi^ 


\ Ms/ 





3 Jd f^a + f^s{y) Ms(y) 

Let No{x, y) be the Neumann function of ^ — 3^ A on Q. We suppose that satisfy the 
ellipticity condition (|2.1[) . It follows from Theorem 12.11 that there is a constant C such that 



\N{x,y) ~ No{x,y)\ < Cfislx - y\-^+\ 
\V,{N{x, y) ~ No{x, y))\ < Cfi.,\x - y\--'+\ 
\Wl{N{x, y) - No{x, y))\<C - y\-'+^ + M^|x - y\-'+^) , 

for all x,y €z D provided that uj is bounded. On the other hand, it is proved in [31 Lemma 
2.2] that there is a constant C such that 

|iVo(x,2/)-3M,r(x-2/)|<CMp, 
|V,(7Vo(x,y) - 3M,r(x -y))\<C (m? + Mp|a; - y\-') , 

\\7l{No{x, y) - 3M,r(x -y))\<C [fi'J' + fj'\x - y\-') , 

for all x,y G D provided that e^/ms is sufficiently small. Therefore, if we put 

Rix,y)^Nix,y)-3fLsrix^y), (3.13) 
we obtain the following lemma. 
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Lemma 3.2 Let R be defined by iS.lSfl . There exists a constant C such that 

\R{x, y)\<C (fj^ + Jis\x - 2/1-1+^) , (3.14) 

|V,i?(a;,j/)| <c(/22+^3/2|^„^|-i^-^|^._y|-2+Aj ^ (3^15) 

\^lR{x, y)\<C (jlV^ + fj^\x - y\-^ + f^M - y\-'^^ + M^k - ^1"'+') • (3.16) 
We introduce some notation following [2]. Let 

n{x):^ / N{x,y)dy, xeD, (3.17) 

JD 

and define a multiplier M by 

M[f]{x)- iian{x)f{x). (3.18) 
We then define two operators M and TZ by 

M[f]{x) := i^ians I ifiy) - f{x))r{x - y) dy 

JD 

+ / {—^ -)V/(y) • VyT{x - y) dy, (3.19) 

nf]{x) Ma / (./(y) - f{x))R{x,y) dy 
Jd 

+ \[ {—^ -)V/(2/) • V,i?(x, y) dy. (3.20) 

3 J_D Ma + Ms Ms 

Then, p.l2p can be rewritten as 

(/-7W)[$] - (A/' + 7^)[$] = onD, (3.21) 

where / is the identity operator. 
For > 0, define 

Then one can show using Holder's inequality that 



||r,[/]|U.(,,)<Ce"||/|U.(z,) (3.22) 



for all p > I . 



We fix A so that A > i. Suppose that e^^/TJ and Using p.22p one can show that 

WmWw^.HD) < C (^eVaMs + l|V/||i,(^) . (3.23) 
One can also show using p.l4p - p.l6p and p.22p that 

\\n[f]\\w^..(D) < Ce^Ml [t^aj^se^ + ^) W^fh^D)- (3.24) 



Therefore, the following lemma holds. 
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Lemma 3.3 Let p > 3. Then there exists a constant C such that i3.23\) and \3.24-^ hold. 

The rest of derivation of p.9p is exactly the same as in [5] . But we briefly describe it for 
the readers' sake. From p.2ip . we get 

oo 

$ = ^ ((/ - M)-\N + n)y {I - (3.25) 

which converges in W^''p{D) for all p > 3. We then obtain 

'^{x)^{I-M)-^['^^^\x) + {M + n){I^M)-^[<^^'^\x)+E{x), xeD, (3.26) 
where the error term E satisfies 



||£'||w-i,P(_D) < Cep.slla{l + (VJQ [e^^latl■s + II II WLp (D) ■ 



(3.27) 



Then (|3.9p follows from p.26p and the error of the approximation is bounded by a constant 
times 



We emphasize that approximation p.9p is valid under the assumption that e^/^ and ^ 
are small, which indicates that the size and absorption coefficient of the anomaly are much 
smaller than the scattering coefficient. 
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